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Recent investigations on the bifurcations in switching circuits have shown that many atypical bifurcations 
can occur in piecewise smooth maps which can not be classified among the generic cases like saddle-node, 
pitchfork or Hopf bifurcations occurring in smooth maps. In this paper we first present experimental results 
to establish the need for the development of a theoretical framework and classification of the bifurcations 
resulting from border collision. We then present a systematic analysis of such bifurcations by deriving a 
normal form — the piecewise linear approximation in the neighborhood of the border. We show that there 
can be eleven qualitatively different types of border collision bifurcations depending on the parameters of the 
normal form, and these are classified under six cases. We present a partitioning of the parameter space of the 
normal form showing the regions where different types of bifurcations occur. This theoretical framework will 
help in explaining bifurcations in all systems which can be represented by two dimensional piecewise smooth 
maps. 


1 Introduction 


The map is given by 


Most studies in bifurcation theory have been done using 
smooth dynamical systems like the Henon map, the Ikeda 
map and the pendulum equation. In the class of non- 
smootli systems, maps with square root singularity have 
been studied extensively @•110 because of their appli¬ 
cation in impact oscillators and other impacting mechan¬ 
ical systems. On the other hand, piecewise smooth maps 
with finite one-sided partial derivatives at the discontinu¬ 
ity have attracted relatively little attention. Though the 
possibility of strange bifurcations like period-2 to period- 
3 or period-2 to 18-piece chaotic attractor have been re¬ 
ported ||, no systematic study has been made to cate¬ 
gorize the possible bifurcations in piecewise smooth maps. 
Such maps were considered to be just a mathematical pos¬ 
sibility as no physical system with these characteristics was 
known. 

However in recent years there has been a discovery 
that a large class of engineering systems, particularly the 
switching circuits used in power electronics, yield piece- 
wise smooth maps under discrete modeling, and border 
collision bifurcations are quite common in such systems 
H [fj . This has provided motivation for the present study 
whose objective is to systematically analyze all different 
kinds of bifurcations that can occur in two dimensional 
piecewise smooth maps. 

We consider a general two-dimensional piecewise 
smooth map g(x, y ; p) which depends on a single parame¬ 
ter p. Let T p , given by x = h{y, p) denote a smooth curve 
that divides the phase plane into two regions Ra and Rb- 


g(x,y;p) 


9i(x,y-,p ) for x,y £ Ra, 
92(x,y;p) for x,y £ R B 


(1) 


It is assumed that the functions g± and g 2 are both con¬ 
tinuous and have continuous derivatives. The map g is 
continuous but its derivative is discontinuous at the line 
T p , called the “border”. It is further assumed that the 
one-sided partial derivatives at the border are finite. We 
study the bifurcations of this system as the parameter p 
is varied. 

If a bifurcation occurs when the fixed point of the map 
is in one of the smooth regions Ra or R B , it is one of 
the generic types, namely, period doubling, saddle-node 
or Hopf bifurcation. But if a fixed point collides with 
the borderline, there is a discontinuous jump in the eigen¬ 
value of the Jacobian matrix. In such a case an eigen¬ 
value may not “cross” the unit circle in a smooth way, 
but rather “jumps” over it as a parameter is varied con¬ 
tinuously. One therefore cannot classify the bifurcations 
arising from such border collisions as those occurring for 
smooth systems where the eigenvalues cross the unit circle 
smoothly. In this paper we develop a new classification for 
border collision bifurcations. 

The paper is organized as follows. In Sec. 2, we illus¬ 
trate the problem with the help of an example of switching 
circuit. In Sec. 3, the normal form is derived. In Sec. 4, we 
analyse the border collision bifurcations occurring in piece- 
wise smooth maps. We present a partitioning of the pa¬ 
rameter space of the normal form exhibiting various kinds 
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of border collision bifurcations. We conclude in Sec. 5. 


2 Examples of border collision bi¬ 
furcations in a power electronic 
circuit 

The subject of power electronics is concerned with high 
efficiency conversion of electric power, from the form avail¬ 
able at the power source, to the form required by the spe¬ 
cific appliance or load. Power electronic technology is in¬ 
creasingly finding application in the home and workplace: 
familiar examples are domestic light dimmers, fluorescent 
lamp ballasts, battery chargers, and switch-mode power 
supplies of all electronic appliances including the personal 
computer. 

In contrast with mainstream electronics, power electron¬ 
ics is characterized by the use of electronic switches which 
operate in “on” or “off” state. Since electrical power sup¬ 
plies can be either dc or ac, there are four basic types of 
power converters: ac-dc, dc-ac, dc-dc and ac-ac. Here we 
will consider one of the simplest but most useful of power 
converters — the dc-dc buck converter — which is used to 
convert a dc input to a dc output at a lower voltage. 


Triangular 

wave 



(a) 



Figure 1: (a) The buck converter with duty cycle con¬ 
trolled by voltage feedback, (b) The three ways the state 
can move from one sampling instant to the next. 

The circuit diagram of the buck converter is shown in 
Fig.[l](a). The controlled switch S (generally realized by a 
MOSFET) opens and closes in succession, thus “chopping” 
the dc input into a square wave that alternates between 



Figure 2: Experimental bifurcation diagram of the buck 
converter. The parameter values are: R = 23.5S2, C = 
5/rF, L = 2.96mH, triangular wave: Vu = 8.43V, Vl = 
3.62V, frequency 12 kHz. Bifurcation parameter Vi n varied 
from 35 to 75 V. 

the input voltage V n and zero. The pulsed waveform is 
then low-pass filtered by a simple LC network, removing 
most of the switching ripple and delivering a relatively 
smooth dc output voltage v to the load resistance R. The 
diode D provides a path for the continuation of the induc¬ 
tor current during the off period. The dc output voltage 
can easily be varied by changing the duty ratio, i.e., the 
fraction of time that the switch is closed in each cycle. 

In practice it is necessary to regulate v against changes 
in the input voltage and the load current. For example, if a 
buck converter is used to convert the standard 5 V dc sup¬ 
ply used in computers to the 3.3 V needed for the Pentium 
CPU chip, it would be necessary to regulate the average 
output voltage at 3.3 V in spite of the varying power de¬ 
mand of the chip. This can be achieved by controlling 
the switch S by voltage feedback as shown in Fig.[T|. In 
this simple proportional controller, a constant reference 
voltage V re f is subtracted from the output voltage and 
the error is amplified with gain A to form a control sig¬ 
nal v C on = A{v — V re f). The switching signal is generated 
by comparing the control signal with a periodic sawtooth 
(ramp) waveform. S turns on whenever v con goes below 
Vramp and a latch allows it to switch off only at the end of 
the ramp cycle. 

Though this circuit or its variants are used in a large 
number of practical applications requiring regulated dc 
power supply, it has been demonstrated [|[ |l(J that the 
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Figure 3: Experimental bifurcation diagram of the buck 
converter. The parameter values are: R = 28.9S2, C = 
5/zF, L = 2.96mH, triangular wave: Vjj = 8.43V, Vl = 
3.62V, frequency 8 kHz. Bifurcation parameter varied 
from 50 to 70 V. 

system can exhibit bifurcations and chaos for a large por¬ 
tion of the parameter space. To investigate the dynamics 
analytically, we obtain a two dimensional Poincare map 
by sampling the inductor current and capacitor voltage at 
the end of each ramp cycle. 

Because of the transcendental form of the equations, the 
map cannot be determined in closed form. In simulation, 
the map has to be obtained numerically. It is however 
possible to infer the form of the map. There are three ways 
in which the system can move from one observation point 
to the next: (a) the control voltage is throughout above the 
ramp waveform and the switch remains off, (b) the cycle 
involves an off period and an on period, (c) the control 
voltage is throughout below the ramp waveform and the 
switch remains on. The three cases are shown in Fig. 0(b). 
These are represented by three different expressions of the 
map. The borderlines are given by the condition where 
the control voltage grazes the top and bottom of the ramp 
waveform. Therefore there are three compartments in the 
phase space, separated by two borderlines, and we have a 
piecewise smooth map. 

We present the experimentally obtained bifurcation dia¬ 
grams for this system for different sets of parameter values. 


An experimental bifurcation diagram is shown in 
Fig.|](a). Here we find two parameter values (shown with 
arrows) for which a periodic orbit directly bifurcates into 
a chaotic orbit. Such bifurcations have been reported ear¬ 
lier in ||, |ll], |l3|] . The slight expansion of the attractor 

at the bifurcation point is due to system noise and can be 
ignored in theoretical studies. In Fig.|](b) we present the 
continuous time plots of v CO n and the triangular wave volt¬ 
age at the bifurcation point shown by the second arrow, 
where a period-3 orbit bifurcates into a 3-piece chaotic or¬ 
bit. It is seen that the v con waveform grazes the top of 
the triangular wave, which means that a border collision 
bifurcation has occurred. 

The distinguishing feature of this chaotic attractor is 
that there is no periodic window over a large range of the 
parameter value. We find from simulation that there are 
no coexisting attractors in this range. We say a chaotic 
attractor is robust if, for its parameter values there ex¬ 
ists a neighborhood in the parameter space with no peri¬ 
odic attractor and the chaotic attractor is unique in that 
neighborhood |Q. The chaotic attractor resulting from 
this border collision is therefore robust. The question is, 
under what condition does robust chaos occur? 



Figure 4: Numerically obtained bifurcation diagram of the 
buck converter. The parameter values are: R = 22f2, C = 
47/rF, L = 20mH, triangular wave: Vu = 8.2V, Vl = 3.8V, 
time period 400^s. 

Another experimental bifurcation diagram for this sys¬ 
tem is shown in Fig|i|(a). The arrow shows a period dou¬ 
bling bifurcation, but the two bifurcated orbits do not di¬ 
verge perpendicularly from the path of the fixed point be¬ 
fore the critical parameter value. This is therefore not a 
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standard pitchfork bifurcation. This kind of bifurcation 
has been reported in |l5|, [l(| also. Fig.|j(b) gives the con¬ 
tinuous time plots of v con and the triangular wave voltage 
just after the bifurcation and shows that the period dou¬ 
bling occurred at a border collision. Again the question 
is, under what condition does this special type of period 
doubling occur? 

It has been reported earlier 0 that this system has 
coexisting attractors for some ranges of parameter values. 
Since multiple attractors cannot be seen in experimental 
bifurcation diagrams, we present a numerically obtained 
bifurcation diagram in Figg showing the evolution of the 
main attractor and a coexisting attractor. It is found that 
the chaotic attractor comes into existence out of nothing 
at a particular parameter value. Under what condition 
can such strange bifurcations occur? 

In the following sections we develop a complete theory 
of bifurcations in piecewise smooth maps, from which the 
answers to the above questions can be derived. 


3 The normal form 


Since the local structure of border collision bifurcations 
depends only on the local properties of the map in the 
neighborhood of the border, we study the border collision 
bifurcations with the help of “normal forms” — the piece- 
wise affine approximations of g in the neighborhood of the 
border. 

Define 

x = x — h(y, p), and y = y. 

This p-dependent change of variables moves the border to 
the y axis. Then the map g(x 1 y; p) can be written 

g(x + h(y; p), y; p) = f(x, y; p), 

and the border is x = 0. Suppose that when p = po the 
map /(i, y; p) has a fixed point Po on the border, that is, 

Po = (0,y o (po)) = /(0,yo(po);Po)- 


Let ei be a tangent vector in the y direction. The vector 
ei maps to a vector e 2 - We assume e 2 is not parallel 
to ei. Define the local coordinates as the following (C.f. 
Fig. g). Choose the point Po as the new origin for ei in 
the y direction and e 2 in the x direction. In these x-y 
coordinates, the fixed point Po is given by (0,0), and the 
border T p is given by x = 0. We define the new parameter 
p = p — po so that po = 0. Choose the scales such that at 
p = 0 an unit vector along y- axis maps to an unit vector 
along ai-axis. The phase space is now divided into the two 
halves L and R and the map f{x,y\p) can be written as 
F(x, y; p). 

We can write the map F(x,y;p) in the side L in the 
matrix form as 


F{x, y; p) 


h(x,y;p) 

f 2 (x,y;p) 


and F(0, 0; 0) 



Linearizing F(x, y; p) in the neighborhood of (0,0;0), we 
have: 


F{x,y;p) 


( Jn Ji 2 W x \ + f v Lx 
\ J 21 J 22 ) \ y ) V Vl v 
o(x, y ; p), for x < 0 


where 

d 

J11 = Ibn —/i(a;,y; 0 ), 

x —> 0 “ , y —>0 OX 

d 

J12 = lim — /i(aj, 2/; 0 ), 

x— » 0 - , y — >0 oy 

d 

J21 = lim irrf2(x,y;0), 

x —> 0 - , y —>0 OX 

d 

J22 = lim ^-f2{x,y;0), 

x^o~ , 0 oy 

d 

v Lx = lim — /i(x,y; 0 ), 

X^ 0 ~,y ^0 oy 

d 

v Ly = lim -^rf2(x,y-,0). 

x —> 0 — , y — >0 op 


+ 


( 2 ) 


Since an unit vector along y- axis maps to an unit vector 
along x-axis, by (g) this particular choice of coordinates 
makes J 12 = 1 and J 22 = 0. Further, we note that Jn 
is the trace (denoted tl) and J 21 is the negative of the 
determinant (denoted —Sl) of the Jacobian matrix. Thus 
© becomes 


F{x,y,p) 



+o(x,y,p), if x < 0, (3) 


Similarly, for side R we obtain 


F(x,y;p ) 



+o(x,y;p), if x > 0. (4) 


where the corresponding quantities in R are defined in a 
similar way. 

Continuity of the map implies 

( v Lx \ _ ( VR X \ 

V v Ly ) ~ V v Ry ) 

We now make another change of variables so that the 
choice of axes is independent of the parameter. The co¬ 
ordinate transformation x = x, and y = y — pv y , and 
p = Jl (y x + v y ) (assuming ( v x + v y ) ^ 0) gives 



G 2 = 


tl 

-5l 


tr 

-5r 


1 

0 

1 

0 


y ) + M ( J ) > for a; < 0, 

y ) + M (0) ’ for * > °’ 
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Figure 5: The transformation of coordinates from the two dimensional piecewise smooth map to the normal form. 


which is the desired 2-D normal form. 

Note that if (■ v x + v y ) = 0, then the fixed point moves 
along the border as p varies. Hence we assume the generic- 
ity condition (v x +v y ) ^ 0 to ensure that a border collision 
occurs at p = 0. 

It is interesting to note that tl and Sl are simply the 
trace and the determinant of the Jacobian matrix of the 
fixed point Po on Ra side of the border P. Let P p denote 
a fixed point of g(x, y ; p) defined on po — e < p < po + e 
for some small e > 0, then P p depends continuously on 
p. Assume that P p is in region Ra when p < po and in 
region Rb when p> po, and that P p is on P when p = po- 
For p < po, the eigenvalues of the Jacobian matrix of the 
fixed point P p are denoted as Ai and A 2 . Since the trace 
and the determinant of the Jacobian is invariant under the 
transformation of coordinates, we can obtain the values of 
tl and Sl as 


t l = lim (Ai + A 2 ) 

P^Po 


S L = lim (AiA 2 ). 

P^Po 

(6) 


The values of tr and Sr can be calculated in a similar way 
for p> po- This property is very important in numerical 
computations. For a border-crossing periodic orbit with 
higher period, we examine the pth (if the period is p ) it¬ 
erate of the map. The matrices in (ra) then correspond to 
the pth iterate rather than the first iterate of the map. 

When Sl and Sr are zero, the system becomes one di¬ 
mensional and the normal form reduces to 


Gi{x;p) 


a x + p for x < 0, 
b x + p for x > 0, 


(7) 


4 Classification of border collision 
bifurcations 

Various combinations of the values of tl, tr. Sl and Sr 
exhibit different kinds of bifurcation behaviors as p is var¬ 
ied through zero. To present a complete picture, we break 
up the four dimensional parameter space into regions with 
the same qualitative bifurcation phenomena. If the pa¬ 
rameter combination is inside a region, then g and G 2 will 
have the same types of bifurcations. If it is on a bound¬ 
ary, then higher order terms are needed to determine the 
bifurcations of g. 

The fixed points of the system in both sides of the 
boundary are given by 


L* 

R* 


A* 


1 - tl + Sl’ 1 


-$LP \ 

- t l + Sl J 


A 4 


1 - Tr + Sr ’ 1 


SrP \ 
- Tr + Sr) 


and the stability of each of them is determined by the 
eigenvalues A] j2 = | (r ± \/t 2 — 45). If the eigenvalues are 
real, the slopes of the corresponding eigenvectors are given 
by — (<5/Ai) and —(S/X 2 ), respectively. Since we consider 
only dissipative systems, we assume |<5 l| < 1 and |<5 r| < 
1. Under this condition there can be four types of fixed 
points. 


1. When S > r 2 /4, both eigenvalues of the Jacobian are 
complex, indicating that the fixed point is spirally 
attracting. If r > 0, it is a clockwise spiral, and if 
t < 0 the spiralling motion is counter-clockwise. 


where a and b are the slopes of the graph at the two sides 
of the border x = 0. 
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2. When S < r 2 /4, both eigenvalues are real. If 2 \/5 < 
t < (1 + J) then the eigenvalues are positive and the 
fixed point is a regular attractor. If —2 \fS > r > 
















— (1 + 5) then the eigenvalues are negative and it is a 
flip attractor. 

3. If t > 1 + <5, then 0 < A 2 < 1 and Ai > 1. The fixed 
point is a regular saddle. 

4. If r < — (1 + 5), then A 2 < — 1 and — 1 < Ai <0. The 
fixed point is a flip saddle. 

If the determinant is negative, there can be only two 
types of fixed points: 

1. For — (1 + £)< t <(1 + <5), one eigenvalue is positive 
and the other negative — which means that the fixed 
point is a flip attractor. 

2. For r > (1 + 5), Ai > 1 and —1 < A 2 < 0, i.e., the 
fixed point is a flip saddle. If r < —(1 + S), then 
A 2 < — 1 and 0 < Ai < 1. The fixed point is again a 
flip saddle. 

When referring to sides L and R, these quantities have 
the appropriate subscripts, i.e., Ail, A 2 l are the eigenval¬ 
ues in side L and Xi R , X 2 R are the eigenvalues in side R. 
As a fixed point collides with the border, its character can 
change from any one of the above types to any other. This 
provides a way of classifying border collision bifurcations. 

It may be noted that in some portions of the parameter 
space there may be no fixed point in half of the phase 
space. For example, the location of L* calculated by the 
above formula may turn out to be in side R. In such cases, 
the dynamics in L is determined by the character of the 
“virtual” fixed point. We denote such virtual fixed points 
by the overbar sign, as L* and R*. If the eigenvalues are 
real, invariant manifolds of these virtual fixed points still 
exist and play an important role in deciding the system 
dynamics. 

It should also be noted that if a certain kind of bi¬ 
furcation occurs when p is increased through zero, the 
same kind of bifurcation would also occur when p is de¬ 
creased through zero if the parameters in L and R are 
interchanged. Therefore, there exists a symmetry in the 
parameter space and in the following discussion it suffices 
to describe the bifurcations in half the parameter space. 
Moreover, we first consider the case of positive determi¬ 
nant, which constitute a large class of physical systems. 
We take up the special features of systems with negative 
determinant at a later stage. 

A special feature of the normal form (|) is that the 
unstable manifolds fold at every intersection with the x- 
axis, and the images of every fold point is a fold point. 
The stable manifolds fold at every intersection with the y- 
axis and the pre-image of every fold point is a fold point. 
The argument is as follows. Forward iterate of points on 
the unstable manifold remain on the same manifold. In 
the normal form, points on the y-axis map to points on 
the x-axis. As an unstable manifold crosses the y-axis, 


one linear map changes to another linear map. Therefore 
the slope of the unstable manifold in the two sides of the 
x-axis cannot be the same unless the parameters of the 
normal form in the two sides of the border are the same 
(implying a smooth map). In case of the stable manifold, 
the same argument applies for the inverse map. Under 
the action of G 2 , the line x = 0 maps to the line y = 0. 
Therefore under the action of G^ 1 , points on the x-axis 
map to points on the y-axis, and hence the stable manifold 
must have different slopes in the two sides of the y-axis. 

We now present the partitioning of the parameter space 
as shown in Fig|| The system behavior in the various re¬ 
gions of the parameter space are taken up in the following 
subsections. 

4.1 Border collision pair bifurcation 

Case 1: If 

t l >(1 + <5l) and t r <(1 + < Sr), (8) 

then there is no fixed point for p < 0 and there are two 
fixed points, one each in L and R, for p > 0. The two 
fixed points are born on the border at p = 0. We call this 
a border collision pair bifurcation. An analogous situation 
occurs if tl < (1 + Sl ) and t r > (1 + S R ) as p is reduced 
through zero. Due to the symmetry of the two cases, we 
consider only the parameter region (|^). There can be three 
types of border collision pair bifurcations depending on the 
character of the orbits for p > 0. 

Case 1(a): If (1 + 6 R ) > t r > —(1 + Sr), then R* is 
stable. Therefore it is like a saddle-node bifurcation, where 
a periodic attractor appears at p = 0. There are two special 
features of this saddle node bifurcation. First, the fixed 
points are born on the border and move away from it as p 
is increased. Second, there is no intermittency associated 
with this bifurcation. 

Case 1(b): If 

t l >( 1+S l ) and r R <-(1 + S R ) and (9) 
SltrXil—Sr\ilX2l+5 r \2l— Sltr+tlSl —S\ — XilSl > 0 

, (10) 

there is a bifurcation from no attractor to a chaotic at¬ 
tractor. The chaotic attractor for p > 0 is robust Q. 
Case 1(c): If tl > (1 + Sl) and tr < — (1 + S R ) and 

SltrXil—Sr\ilX2l+Sr\2l—Sltr+tlSl—S\— \2lSl < 0 

then there is an unstable chaotic orbit for p> 0. 

For (|)|). L* is a regular saddle and R* is a flip saddle. 
Let U l and Sl be the unstable and stable manifolds of L* 
and Ulj and S R be the unstable and stable manifolds of 
R* , respectively. As shown earlier, Ul and U/j experience 
folds along the x-axis, and all images of fold points are fold 
points. Sl and S R fold along the y-axis, and all pre-images 
of fold points are fold points. 
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Figure 6 : The partitioning of the parameter space into regions with the same qualitative bifurcation phenomena. The 
numbering of the regions are the cases as discussed in the text. The regions shown in primed numbers have the same 
bifurcation behavior as the unprimed ones when /r is varied in the opposite direction. 


For condition (§)> Ail > A 2 l > 0 and 0 > Ai/{ > A 2 r. 
The stable eigenvector at R* has a slope mi = (— Sr/Xir) 
and the unstable eigenvector has a slope m2 = (—Sr/X2r). 
Since points on an eigenvector map to points on the same 
eigenvector and since points on the y-axis map to the x- 
axis, we conclude that points of to the left of y-axis 
map to points above x-axis. From this we find that U/j has 
an angle m 3 = ( 5 lX 2 r)/( 6 r — tlX 2 r) after the first fold. 
Under condition (|9j) we have mi > m 2 > 0 and m 3 < 0. 
Therefore there must be a transverse homoclinic intersec¬ 
tion in R. This implies an infinity of homoclinic intersec¬ 
tions and the existence of a chaotic orbit. 

We now investigate the stability of this orbit. The basin 
boundary is formed by Sl- Sl folds at the y-axis and 
intersects the x-axis at point C. The portion of Ul to 
the left of L* goes to infinity and the portion to the right 
of L* leads to the chaotic orbit. Ul meets the x-axis at 
point D , and then undergoes repeated foldings leading to 
an intricately folded compact structure as shown in Fig.^j. 

The unstable eigenvector at L* has a negative slope 
given by (— Sl/Xil)- Therefore it must have a hetero¬ 
clinic intersection with S#. Since both Ul and U# have 
transverse intersections with Sr, by the Lambda Lemma 


JT| we conclude that for each point gonU/j and for each 
e-neighborhood N e (q), there exist points of Ul in N e (q). 
Since Ul comes arbitrarily close to U#, the attractor must 
span Ul in one side of the heteroclinic point. 

Since all initial conditions in L converge on Ul and all 
initial conditions in R converge on U#, and since there 
are points of Ul in every neighborhood of \ 1 r, we con¬ 
clude that the attractor is unique. This chaotic attractor 
cannot be destroyed by small changes in the parameters. 
Since small changes in the parameters can only cause small 
changes in the Lyapunov exponents, where the chaotic at¬ 
tractor is stable, it is also robust. 

It is clear from this geometrical structure that no point 
of the attractor can be to the right of point D. If D lies 
towards the left of C , the chaotic orbit is stable. If D falls 
outside the basin of attraction, it is an unstable chaotic 
orbit or chaotic saddle. From this, the condition (|n]) of 
stability of the chaotic attractor is obtained. If Sl =Sr = 5 
this condition reduces to trXil — Ail A 2 L + tl — tr — 5 > 0 . 
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Figure 7: The stable and unstable manifolds of L* for 
tr = 1.7, 5l = 0.5, tr = —1.7, Sr = 0.5. R* is marked by 
the small cross inside the attractor. 

4.2 Border crossing bifurcations 

In all regions of the parameter space except (|j), a fixed 
point crosses the border as p is varied through zero. The 
resulting bifurcations are called border crossing bifurca¬ 
tions. In the following discussions we consider the bifur¬ 
cations as p varies from a negative value to a positive 
value. 

Case 2: Linear attractor to flip saddle. This occurs 

if 2 \ff L < t l < (1 + S L ) and t r < -(1 + Sr). 

There is a bifurcation from a period-1 attractor to a 
chaotic attractor as /a is increased through zero. This 
chaotic attractor is robust. 

For p < 0, L* is a linear attractor while R* is a flip sad¬ 
dle. All initial conditions in L converge on to L*, while 
initial conditions in R converge on to Ur. Since Ur must 
have a heteroclinic intersection with one of the stable man¬ 
ifolds of L, all initial conditions in R also converge on to 
L*. 

For p > 0, R* is a flip saddle. As shown in the dis¬ 
cussion for Case 1(b), there is a homoclinic intersection in 
R implying the existence of a chaotic orbit. As L* is in 
R, its stable manifolds point towards R. Since there is an 
intersection of Sr with the invariant manifold associated 
with Air, all initial conditions converge on Ur, making 
the chaotic attractor unique. 

Case 3: There is a unique period-1 attractor for both 
positive and negative values of p in the following cases. At 
border collision, only the path of the fixed point changes. 


Regular attractor to spiral attractor: This occurs 

if 2\/~Sl < t l < (1 + Sl) and - 2^/Sr < tr < 2^/Sr. 

For p < 0, all initial conditions in R are attracted to 
R* which is in L. All initial conditions in L converge on 
to L*. Therefore the fixed point is the unique attractor. 
For p > 0, all initial conditions in L move linearly towards 
L* which is in i?, and all points in R spiral towards R*. 
Therefore R* is the unique attractor. 

Spiral attractor to spiral attractor having the same sense 
of rotation: This occurs 

if 0 <t l < 2^[5l and 0 < tr < 2^/5r 

or - 2 \/Sl <t l < 0 and - 2 \/Sr < t r < 0. 

If the spiralling orbits in L and R have the same sense, 
there is an overall spiralling orbit converging on the fixed 
point. Therefore there is an unique period-1 attractor for 
both p < 0 and p > 0. 

Regular attractor to regular attractor: 

2\J~Sl < t l < (1 + Sl) and 2^/Sr < tr < (1 + Sr), 

Flip attractor to flip attractor: 

-2> tl > -(1+5r) and -2 <J8r > tr > — (1+<5r), 

Regular attractor to flip attractor: 

2\[Sl < tl < (1 + Sl) and - 2> tr > -(1 + Sr). 

In the above three cases, for p < 0, initial conditions in 
R move linearly to R* . Since there must be a heteroclinic 
intersection of the stable manifolds, all initial conditions 
converge on L* . The situation for p > 0 is similar. 

Case 4: In the following cases there can be bifurcation 
from multiple attractors to multiple attractors. There are 
general mechanisms for the occurrence of coexisting at¬ 
tractors. 

Spiral attractor to spiral attractor with opposite sense of 
rotation: This occurs 

if 0 < tr < 2\Ar and — 2^/Sr < tr < 0 

or - 2 \JJl < tr < 0 and 0 < tr < 2*JFr. 

Spiral attractor to flip attractor: This occurs 

if —2\/<5 r < tr < 2^/Sl and -2 ^/Sr > t r > -(1+<5r). 

There can be multiple attractors on both sides of p, one 
of which is a fixed point. 

Case 5: In the parameter space region 

tr < -(1 + Sr) and tr < 0, 





initial conditions in L move to R and vice versa. Therefore 
the dynamics is governed by the stability of the second 
iterate with one point in L and the other in R. 

The eigenvalues of the second iterate are 

\ ( T L Tr — Sr — S L ± 

\Jt l 2 t r 2 - 2 t l tr6r - 2 t l S l t r + S R 2 - 2 SrSl + Sl 2 ^ 

From this, the condition of stability of the period-2 orbit 
is obtained as 


conditions diverge away from it along the unstable eigen¬ 
vector. Therefore there can be no attractor for p > 0. 

Case 5(d): If 

t l <-(1 + S l ) and t r < -(1 + Sr), 

there is no attractor for both positive and negative values 
of p since all the fixed points of the first and second iterate 
are unstable. 

Case 6 -.Spiral attractor to flip saddle: This occurs 
if 0 <t l < 2 and t r < -(1 + Sr). 


1 - t l t r + S l + Sr + S l Sr > 0 for Ai<+l(ll) 
1 + t l t r - S L - S R + S L S R > 0 for A 2 >-1(12) 

There are three sub-cases: 

Case 5(a): If 

tr < -(1 + Sr) and t l < -2 y/Sf 
and 1 - t l t r + S l + Sr + S l Sr > 0, 

then there is a unique period -1 attractor for /i < 0 and a 
unique period -2 attractor for p > 0 . 

For p < 0, L* is a flip attractor and R* is a flip sad¬ 
dle. All initial conditions in L converge on L* and all 
initial conditions in R go to L in the first iteration and 
then converge on to L*. For p > 0, the condition (0 
ensures the stability of the period-2 orbit. The existence 
of heteroclinic intersection makes the attractor unique. 

This is like a period doubling bifurcation occurring on 
the borderline. In contrast with standard period doubling 
bifurcation, the distinctive feature of the border collision 
period doubling is that as p is varied through zero, the 
bifurcated orbit does not emerge orthogonally from the 
orbit before the bifurcation. 

Case 5(b): If 

tr < -(1 + Sr) and -2*/5f<T L <0 
and 1 - t l t r + S L + Sr + 5 l 5r > 0, 

then for p < 0 there can be multiple attractors one of 
which is a period-1 fixed point. For p > 0, the period-2 
orbit involving both L and R is stable. Therefore there is 
an unique period -2 attractor. 

Case 5(c): If 

tr < -(1 + Sr) and - (1 + S L ) < t l 
and 1 - t l t r + S L + Sr + S l Sr < 0, 

then there is a period-1 attractor for p < 0. For 
— (1 4 - Sl) < tl < — 2 \/5l, the eigenvalues of L* are real 
and coexisting attractors cannot occur. Otherwise mul¬ 
tiple attractors may exist. For p > 0, since ( 0 ) is not 
satisfied, it implies that the fixed point of the twice iter¬ 
ated map is unstable. The eigenvalues are real and initial 


For p < 0, there can be multiple attractors one of which is 
a period-1 fixed point. The asymptotic behavior for p> 0 
may be a periodic attractor (of periodicity greater than 
unity), or chaotic attractor. As Tj J is increased, periodic 
windows of successively higher periodicities (2,3,4,...) oc¬ 
cur, and there are windows of chaos between two such peri¬ 
odic windows. The period-n attractor comes into existence 
through a border collision pair bifurcation in the nth iter¬ 
ate and goes out of existence when the period-n fixed point 
becomes unstable. From the stability boundary of 

period -2 attractor is given by 1 +tlTr—Sl~Sr+SlSr = 0 . 
For higher iterates such analytical expressions for the 
boundary of periodic windows become involved and are 
not presented here. There is no mechanism to prevent the 
occurrence of multiple attractors. 

This gives a complete description of the bifurcations 
that can occur at various regions of the parameter space of 
the normal form (|^) . Representative bifurcation diagrams 
of the cases (where attractors exist) are shown in Fig.|j. 

4.3 The case of negative determinant 

If the determinant is negative, one has to find out which 
type of fixed point changes to which type at it moves across 
the border. Depending on the type of the fixed point at the 
two sides of the border, the bifurcations will be of the same 
kind as discussed in the previous section. For example, if 
Sl,Sr < 0 then the eigenvalues are real for all values of tl 
and tr. Therefore there can be no coexisting attractors 
anywhere in the parameter space. The region of stability 
of period-2 attractor, given by conditions (12) and (13), 
is much larger. Moreover there is a region of parameter 
space where a border collision pair bifurcation results in 
the creation of a period-2 attractor since condition (13) 
is satisfied. The partitioning of the parameter space for 
negative determinants is given in Fig.9. 

There is however a difference in the equation for the 
boundary crisis in border collision pair bifurcation. For 
— 1 < Sr < 0, we have 1 > Aik > 0, A 2 R < —1, and R* is 
located above the x-axis. A positive value of Aik implies 
that U l converges on Uk from one side. If 

Ai y‘ > Az y \ as) 

t l 1 Sl tr-I-Sr 
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Figure 8: Representative bifurcation diagrams of the normal form when fj, is varied from a negative value to a positive 
value. For the cases where multiple attractors can exist, only one of many possibilities is shown. Case 1(a): No 
attractor to period-1 attractor; Case 1(b): No attractor to chaos; Case 2: Period-1 to chaos; Case 3: Period-1 to 
period-1; Case ): Period-1 + period-3 coexisting attractors to period-1 + period-4 coexisting attractors; Case 5(a): 
Period-1 to period-2; Case 5(b): Period-1 + period-11 coexsting attractors to period-2 attractor; Case 6: Period-1 to 
coexisting period-5 + chaotic attractors. 


then the intersection of \5 l with the x-axis remains the 
rightmost point of the attractor and (]l(j|) still gives the 
parameter range for boundary crisis. But if ( |l3| ) is not 
satisfied, the intersection of with the x-axis becomes 
the rightmost point of the attractor, and the condition of 
existence of the chaotic attractor changes to 

Mr — 1 Sl (jl — Sl - Ml) 

tr-I-Sr (tl — 1 — Sl) {SrMl — 8ltr) 

For Sl <0 and Sr < 0, L* is below the x-axis and the 
same logic as above applies. But if <5^ < 0 and Sr > 0, the 
stable manifold of R* has a negative eigenvalue and hence 
U l does not approach Ufj from one side. Therefore, if ( |l3| ) 
is not satisfied, there is no analytic condition for boundary 
crisis — it has to be determined numerically. 

5 Conclusions 

In this paper we have investigated the various types of 
border collision bifurcations that can occur in piecewise 
smooth maps by deriving a piecewise affine approximation 
of the map in the neighborhood of the border. We have 
shown that there can be basically eleven different types of 


border collision bifurcations, classified under six “cases”. 
We have presented a partitioning of the parameter space 
into regions where qualitatively different bifurcations oc¬ 
cur. 

This body of knowledge helps us in explaining the bi¬ 
furcations observed in experimental and numerical inves¬ 
tigations of switching circuits, some of which have been 
presented in Sec.2. For example, the experimental bifur¬ 
cations of the type seen in Fig.|| can occur in Case 2 and a 
part of Case 6. A period doubling bifurcation of the type 
shown in Fig.|^ can occur in the second iterate of the map 
if the parameters fall under Cases 5(a), 5(b) and a part of 
Case 6 (coexisting attractors can not be observed in ex¬ 
perimental bifurcation diagrams). The sudden appearance 
of a chaotic attractor as in Fig.^ can occur in border col¬ 
lision pair bifurcation and can be categorized under Case 
1 (b). Note that this bifurcation occurs in the third iterate 
while the period-1 attractor is present, and therefore the 
resulting chaotic attractor is not robust. 

The theoretical problem dealt in this paper was posed 
by the recent investigations in switching electrical circuits. 
But we believe that such atypical bifurcations will be ob¬ 
served in other nonsmooth physical systems also and the 
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Figure 9: Schematic diagram of the parameter space par¬ 
titioning for —1 < 8 l < 0 and — 1 < Sr < 0 into regions 
with the same qualitative bifurcation phenomena. (1) No 
fixed point to period-1; (2) No fixed point to period-2; (3) 
No fixed point to chaos; (4) No fixed point to unstable 
chaotic orbit, no attractor; (5) Period-1 to period-2; ( 6 ) 
Period-1 to chaos; (7) Period-1 to period-1; ( 8 ) Period-1 
to no attractor; (9) No attractor to no attractor. The re¬ 
gions shown in primed numbers have the same bifurcation 
behavior as the unprimed ones when p is varied in the 
opposite direction. 


theory developed in this paper will help in understanding 
the nonlinear phenomena and bifurcations in such systems. 
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